Transition temperature of the homogeneous and dilute Bose gas in D-dimensions 
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The phase transition temperature of the homogeneous and dilute Bose gas in D-dimensions (2 < 
D < 3) is calculated by a mean field-based statistical method. The shift of the phase transition 
temperature is written up to the leading order as AT C /T® = cry a , where 7 = n 1,/3 a. We derived 
Huang's result of the phase transition temperature in the generalized dimensions. We show that 
c(D) is positive and a(D) = 2(15/2 — l) 2 in the short-wavelength range. The origin of the difference 
between a — 1/2 and a — 1 at D=3 is discussed. The T c at D = 2 is calculated in the same scheme. 
The result is compared with Fisher and Hohenberg's KT temperature. 
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I. INTRODUCTION 

The homogeneous and dilute(HD) Bose gas is a funda- 
mental topic of physics to which a lot of theoretical ef- 
forts have been devoted in recent years. The shift of the 
transition temperature(T c ) remains an interesting prob- 
lem to this day. At the same time, the study of T c in 
two-dimensions(2D) has been a crucial topic because it 
is directly connected with the Kosterlitz-Thouless(KT) 
transition. 

In a homogeneous system, the density in D-dimensions 
is expressed as n(D) = n(l) D . The diluteness is ex- 
pressed by a dimensionless gas parameter 7 = n(l)a = 
n(3) 1 / 3 a, where n(l) is the one-dimensional number den- 
sity and a is the s-wave scattering length. In a repulsive 
and dilute gas, 7 is positive and much less than 1. The or- 
der of diluteness can be defined as | I117I, and much larger 
than 1. The argument in the paper is valid only for an 
extremely weakly interacting Bose system, and therefore 
cannot be applied to a strongly interacting system such 
as real liquid 4 He. 

The transition temperature of the non-interacting ho- 
mogeneous Bose gas in D-dimensions is [l[ 



T° C {D) = 



2-Kh 2 n(2) 
mk B C(§)ts' 



(1) 



where 2 < D < 3 and C is the Riemann-zeta function. 
The transition temperature of the HD Bose system is 
written as a series form of 7. The leading orders in 3D 
and 2D are 



T c (3) = r c °(3)(l + C 7 Q + ■■■), 

_/ 

T c (2) - T c °(3)(0- 



hi7| Q ' 



(2) 
(3) 



The difficulty in finding the c comes from a break- 
down of perturbation theory at the second order tran- 
sition temperature. Comparing with Eq. Ip]). the inter- 
acting system is equivalent to using the effective mass 



m* = m/(l + cj a ) in the ideal gas expression for T°(3). 
The positive c means that the effective mass of the quasi- 
particle is reduced and vice versa for the negative c. 

Two choices of a for the leading order have been deter- 
mined with a repulsive hard-sphere interaction a = 1/2, 
suggested by Toyota and Huang from mean field-based 
microscopic theories 

ill. 

On the other hand, a = 1 
was suggested by others 0, H, ■ Toyota and Huang 
reported the coefficient c as 



8V2^ 



3.524, 



(4) 



but the others claimed that c is on the order of 1. Nowa- 
days, a = 1 and positive c has been widely accepted since 
numerical simulations @, Q . 

Fisher and Hohenberg(FH)'s paper has suggested a 
KT transition temperature in 2D for the HD Bose gas 
from the Bogoliubov theory under the condition that 
lnln7~ 2 3> 1. It has a double logarithmic form as [l(| 



2ttH 2 ti(2) 
mkB lnln7~ 



(5) 
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At a glance, FH's expression looks different from the Eq. 
©■ 

In this paper, we will calculate the a(D) and c(D) of 
the HD Bose gas between 2 and 3 dimensions with a 
semi-classical method. We will explain why Toyota and 
Huang's leading order a = 1/2 is different from others 
where a = 1. The KT transition temperature in 2D will 
be obtained under the same scheme and compared with 
Hohcnberg and Fisher's, too. 



II. EQUATION OF STATES IN D-DIMENSIONS 

For a momentum-independent interaction, the 
mean field contribution to the self-energy E(-D) is 
7i(D)[/(i9)[llJ, where U(D) is the D-dimensional cou- 
pling constant. Then, the dispersion relation of the 
weakly interacting Bose gas in D-dimensions from the 



2 



mean field theory is 



e = e + E(D), 



(6) 



where £o = p 2 /2m. It corresponds to a short-wavelength 
range in the Bogoliubov energy spectrum [ll[ , and it help 
us to calculate T c without any breakdown problem. In 
principle, this scheme draws from Hartree-Fock theory, 
which assumes that the atoms behave as non-interacting 
bosons in a self-consistent mean field [H, [T3| ■ We may 
choose the self-energy as the exact form 2nU instead of 
nU because the two mean field contributions to the self- 
energy from Hartree and Fock are equal. However, for 
convenience we will use Eq. ^ because the factor of 2 
is out of focus here. 

The U (D) is obtained from a D-dimensional scattering 
theory. Under the hard-sphere interaction of the diame- 
ter a, it is well-known in 3D as 



£7(3) 



Anh 2 a 



(7) 



In 2D. on the other hand, it has a logarithmic form [14j, 

M 



[7(2) = 



2vrft 2 



1 



4irh 2 



1 



m ln(l/fea) m | In -y 2 1 



(8) 



However, we had better use the D-dimensional interpo- 
lated form of U{D) between two and three dimensions 
for general purposes [lj| • 



U(D) = 



^h 2 a D ~ 2 

r(#-i)m 



(9) 



The average occupation number of the momentum- 
independent potential is written from Eq. © as 



P z -l e P(e (p)+S) - 1 

1 



z e - 1 e^o( P ) _ i' 



(10) 



_ -Pi* 



where (3 — 1/ksT. The z is the fugacity given by . 
and 1 below T c . The fj, is the chemical potential, which 
is below the T c and negative above the T c . The z e is 
the effective fugacity given by 



a -/3(|/x|+£) 



(11) 



Note that < z e < 1. The effective chemical potential 
fi e is written in the same way : [i e = [i — S. Therefore, 
the interacting system has the same formula as the ideal 
system except for the effective chemical potential and the 
effective fugacity. 

Let us define the ratio between the effective fugacity 
and the fugacity as 



V 



Ze 
Z 



-/3S 



(12) 



Note that < rj c < 1 and 1 for the ideal system. As 
T — > T c , then z e — > e - ^ 5 " = r\ c . When the interaction 
is very weak, the phase transition is described by the 
variable : r\ c ~ 1 — /3 C S. Therefore, it is clear that even 
for a momentum-independent potential, it shifts the T c . 

From Eq. (fT0|) we have the D-dimensional equation of 
states with z e as [13 



n(D)X D = go(z e ). 



(13) 



A = yj2nh 2 /raksT is the thermal wavelength, and 
g s (z) = YaZi z /I s is the Bose function. In this way we 
can write the transition temperature in D-dimensions, 
2 < D < 3, of the interacting system as 



Tc(D) 



2irh 2 n{2) 
mk B go(ri c )i 



(14) 



From Eq. (|13p the relation of the transition temperature 
between ideal and interacting systems is written as 



(T c ) 2 g»(l) = T c ign( Vc ). 



(15) 



As D approaches 2, the above quantity approaches a fi- 
nite value, 2irh 2 n(2)/mkB, from Eq. (P). It is about 
5.92K for liquid 4 He. 



III. T c BETWEEN 2 AND 3 DIMENSIONS 

From Eq. (fl~5|) . it is necessary to expand the go(rj c ) 
in an analytic form to find the transition temperature. 
Since /3 C S, the interaction term at T c , is much less than 
1, we use the expansion of the Bose functions for D > 2 



i^-'- - , ( f ) +r (i - § j m^-'+oim- 

(16) 



Substituting Eq. (pU into Eq. (|T5|1 . 

D D C(~) 

T r ~ = (T c °)^^^ 



(T°) 



0\ ~2 



C(f) 



■(/3 C S) 



Up to the linear term, we have 



From Eqs. (J) and ^ 

/3 c °(D)Ep) = 



2 r(i-£), n n 



D C(f) 



(17) 



(18) 



mC(f ) D n{l) D ~ 2 U{D) 



2irh 2 



2^-\ 



-\f 1 D\d 

2 / _,D-2 



r(f-i) 



(19) 
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FIG. 1: c(D) as a function of dimensions. 



Note that /?°(3)£(3) = 2( (§) 3 7. The leading order a 
of the shift of the chemical potential E(3)/fes is 1, but it 
definitely cannot be the shift of T c (3). 

Substituting Eq. (flU)) into Eq. |[TB|). we obtain the 
D-dimensional T r as 



T C (D) = T c °(3) 



1 - 



2# 7 r(^- 1 ) 2 r(l-f) 

^C(f)*r(f-i)-- 1 ' 



v 2 (#-l) 2 



(20) 

Therefore, we have the two coefficients of the T c shift in 
D-dimensions as 



c(D) 



2f 7 r(^- 1 ) 2 r(f-f) 

^C(f)*r(f-i)^ 



and 



(21) 



(22) 



where 2 < D < 3. Note that a = in 2D as the Eq. ©. 

There has been a long dispute about the sign of the c. 
From Eq. (J2TJ> we see that it is decided by -T(l - ~ ). 
Therefore, the c is positive at D > 2. It is plotted in 
FIG. 1. However, note that the minimum of c(D) is 
not at D = 2 but at D = 2.17. If D = 3, we obtain 
c(3) = 4\/2^/3e(3/2) 2 / 3 = 1.762 and a = 1/2. Note 
that r(— 5) = -2y/Hr. The c(3) is exactly half of Toyota 
and Huang's microscopic results in Eq. |4]) in that they 
chose 2nU instead of nU as the self-energy. The positive 
c means that the BEC temperature, T c , is different from 
the lambda transition temperature, T\, and there possi- 
bly exists a pseudogap-like region between T c and T\ in 
superfluidity. 



IV. T c IN TWO-DIMENSIONS 

For an ideal Bose system, there is no condensation be- 
cause a long-wavelength phase fluctuation destroys the 
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FIG. 2: The relative transition temperature in 2D compared 
with that of the 3D ideal gas. The solid line is our T c (2) and 
the dashed line is T C FH (2) by Fisher and Hohenberg. 



long-range order. On the other hand, in an interact- 
ing system, the interparticle interactions drive a phase 
transition to a condensate state and then, there exists a 
nonzero condensate density above a specific temperature 
T C KT (2). 

We can also obtain the transition temperature in 2D 
from Eq. ((HJ as 



T c (2) = 



2Trh 2 n(2) 
mk B gi(r)c)' 



(23) 



The convergence of gi(rj c ) is extremely slow. Therefore, 
we need to use another method for the limiting value. 
The Euler's constant C = 0.5772... is defined as 



lim 

S->1 



C(s) 



1 



1 



C. 



(24) 



It is effective up to the first order even when 
limx—,1 gi{x). Then, the limiting behavior of the gi{rj c ) 
when r] c is close to 1 is 



9i(Vc 



1 



1 



\-T) c /3o(2)E(2) 
Substituting Eq. ([25]) into Eq. (|23|) . 

2-nh 2 n{2) 2 U{2) 



T c {2) 



mk 2 B T c (2) 



(25) 



(26) 



Substituting Eq. © into Eq. $So§, we obtain the 2D 
transition temperature as 



T c {2) 



2nh 2 



i(2) 



mk B y/\ l n7 | 



(27) 



The leading orders in Eq. (J3j> are obtained as a' = |, 

and d = C (§)' = 1-897. T c (2) is shown in FIG. 2, and 
compared with Tf H (3) and T c °(3) in the range of interest. 



4 



FH follows the approximation for Tf H (2) under the 
ansatz : /i ~ n/\ ln(a 2 /i)| [101 ]. Then, the transition tem- 
perature in 2D is proportional to the following : 



T C FH (2) 



n(2) 



n(2) 



n(2) 



|lnfea| ln(^a 2 ) In In 7 



(28) 



We think the ansatz is not necessary, because FH's tran- 
sition temperature is nothing but Tf H (2) = E(2)/2fcs 
or just the chemical potential in 2D. Note that the factor 
of half comes from choosing the self-energy nil instead 
of 2nU . This happens because the clean form of U(2) in 
Eq. ([8|) was not known at that time. 



SUMMARY 



Huang's results are just the special case of our calculation 
at D = 3. It was obtained not because the mean field 
theory does not work, but because the long-wavelength 
limit is not included. Therefore, for an extremely low 
density and at the short- wavelength range, we may ac- 
cept Toyota and Huang's results of a = 1/2 [3, 0|. 

The KT transition temperature in 2D in the same con- 
dition is obtained. The two coefficients of the leading 
order of T c (2) (a' and c') are obtained, too. We show 
that Fisher and Hohenberg's KT transition temperature 
is just the chemical potential in 2D. The leading orders 
of the T c in 3D and 2D are obtained as a = a' = 1/2 
in the short-wavelength range, but the reliable results 
of a = a' = 1 are expected if we include the long- 
wavelength limit. 



We derived the transition temperature of the HD Bose 
gas between 2 and 3 dimensions in the short-wavelength 
range using a semi-classical method. The two coefficients 
of the leading order oiT c {D) (a and c), where 2 < D < 3, 
are obtained as a function of the dimensions. Toyota and 
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